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■ Abstract 

G\ '. 

, We argue that the energy dependence of the average position of 

CN ' the shower maximum, X max , and of the number of produced muons, 

N^, can be explained by a change, around 10 17 eV, in the energy de- 
l/-) , pendence of the inelasticity K, which decreases with the energy above 

the string percolation threshold. 

a: 

PS. 

D . The chemical composition of cosmic rays at very high energy, say above 

10 17 eV, is a matter of controversy. The reason being that direct measurement 
.> ; is only possible up to ~ 10 15 eV [1], and a tendency for the effective mass 

^ ■ number A to increase with the primary energy E is observed. Data at higher 

energy, in particular from HiRes experiments [2,3], show that consistency 
between existing theoretical models and observation requires the effective 
mass number A to decrease with energy. This may be due to the presence of 
extra-galactic nucleons [4]. 

In order to make our argument simple, we shall use the original Heitler 
idea [5]: the location of the shower maximum is, on the average, related to 
logE,X max ~ log-E, and the number of charged pions or muons is propor- 
tional to E, N 7T + 7T - ~ ~ i? /3 ,/5<l. If one has a nucleus with A nucleons 
colliding in an hadronic collision we shall write, 
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X max ~X \og(E/A), (1) 

and 

N, ~ A(E/Af, (2) 

or 

logJV M ~(l-/3) log A + piogE. (3) 
We further have, for the elongation length, 



dX rnax — 

and for the logiV M dependence on E, 



dlogA 
~ dlogE 



(4) 



"rflog - ^ = ^ ^rflogE ' (5) 

As experimentally [2,3], above 10 17 eV dX max / dlogE is larger and dlogN^/ 
dlogE is slighty smaller, in comparison with lower energies, the conclusion 
is: 

dlogA/ dlogE < , (6) 

i.e., for E>10 17 eV the average mass number A decreases with energy. 

Another aspect to be taken into account is the X max distribution: it is 
larger for lighter nuclei (in particular in the case of protons). This is con- 
firmed by Monte Carlo simulations [2], and is seen in data: as the energy 
increases the X max distribution becomes wider (more proton-like), thus con- 
firming HiRes results. For a detailed and critical discussion see [6]. 

However, as pointed out in [7], the effect of changing A is equivalent to 
the effect of changing the average inelasticity K. In Fig.l we schematically 
present the situation after the first (hadronic) collision. The quantity (1 — K) 
represents the fraction of energy concentrated in the fastest particle, while 
the inelasticity K is the fraction of energy distributed among produced pions. 
This means that, in the spirit of Heitler model, we shall have, for X maX) 

X max ~ X x + X log[(l -K)^-\, (7) 

instead of (1), where X\ is the average depth of the first collision, Xq is the 
radiation length and E a low energy threshold. For n + , n~ or /x production, 
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(1-k) 



<n> 



Figure 1: Schematic representation of the first hadronic collision. Most of 
the energy fraction, 1 — K, is taken by the leading particle. Most of the 
particles, •n ± 1 come from the fraction of energy left, K. 



N w+W - ~ = CKE 
instead of (2). We further have, 



dX r , 



dhgE 



Xn 



rflog(l — K) 
dlogE 



instead of (4), and 



dlogN„ dlogK 



d log E d log E 
instead of (5). The condition (6) becomes now 

dlogK 



+ 1 



dlogE 



< 



(8) 



(9) 



(10) 



i.e., the inelasticity K has to decrease with the energy. 

Essentially, all the existing high energy strong interaction models based on 
QCD, and QCD evolution, predict an increase with energy - not a decrease - 
of the inelasticity K [8]. The same is true for the hadronic generators Sibyll [9] 
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Figure 2: Energy dependence of the average inelasticity K. Conventional 
QCD model results are from [14]: dotted curve is Sibyll simulation; dashed- 
dotted curve is QGS jet simulations; dashed is a model of [14]; full line is our 
percolation model. 

and QGS jet [10], used in cosmic ray cascade analysis. This happens because 
evolution in the energy implies transfer of energy from valence partons or 
strings, or from bare Pomeron diagram, to sea partons or strings, or to multi- 
Pomeron contributions. 

However, in models with percolation of partons or strings, one expects 
the inelasticity K, above the percolation threshold, to decrease with the 
energy [11]. In the framework of the Dual String Model [13] - but we believe 
the argument is more general - what happens at low energy is the transfer 
of energy from the valence strings to sea strings (and K increases), while 
at higher energy the strings start to overlap in the impact parameter and 
a cumulative effect occurs: the length in rapidity of fused strings is larger 
[12,11]. At some stage, close to percolation threshold, the percolating strings 
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Figure 3: a) Average depth of the shower maximum, X max , b) the num- 
ber density of muons, p M (600)m~ 2 . Data are from [2,3]. Our curves are 
from (7) and (8) with X x = 60gcm" 2 , X = 60gcm" 2 , E = 10 7 eV, 
C = lCT^eV^m- 2 . It is assumed that p^(600m) ~ N^. 

take over the valence strings, and from then on K decreases with the energy. 
Percolation is, in fact, a mechanism for generating fast leading particles. 

In Fig. 2 we show the energy dependence of K in the case of the string 
percolation model [11], in comparison with K determined from QCD inspired 
models, without percolation (see, for instance, [14]). 

In Fig. 3, making use of the values of K for the percolation model of Fig.2, 
and of Eqs. (7) and (8), we show our results for the average depth of the 
shower maximum, X max , and for the logarithm of the number of muons, 
in comparison with HiRes and HiRes-Mia data. As a consequence of the 
decrease of K, and of Eqs. (9) and (10), there is at £>10 17 eV a fast increase 
of X max and a slow increase of logiV M , as seen in data. 

It should be mentioned that low inelasticity events generate large values of 



5 




Figure 4: The relative depth X max — X± as a function of event elasticity 
(1 — K) from [15] by using Sibyll and QGS jet simulations. The curve, from 
(7), for E = 10 19 eV, shows the positive X max — Xi, (1 — K) correlation. 

X max and larger X max fluctuations. This is seen in Fig.4, for X = X max —Xi, 
as a function of event elasticity, taken from [15], by using Sibyll and QGS jet 
simulations incorporated in the CORSIKA program [16]. One further notes 
that X, at fixed energy, has an approximately log behaviour, shown in the 
figure, similar to the behaviour of X as a function of the average inelasticity 
K (Kq-(7)). 

In conclusion, the effect of a decrease of the inelasticity K, for energies 
of the order of 10 17 eV, is similar to the effect of a decrease of the average 
mass number A, and can explain data: Fig.3. We do not claim that there is 
not an increase in the fraction of protons in the primary, what we claim is 
that string percolation, requiring a decrease of K, can qualitatively, at least, 
describe the data around and above E ~ 10 17 eV. 

The reasonable question to ask now is the following: is there a difference 
between the change in A and the change in K7 In our simplified approach, 
mostly based on energy conservation arguments, qualitative differences are 
not seen. However, in the case of percolation the fast particles are not just 
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protons, they may be fast heavy quark (s,c,b,t) mesons with very different 
signatures from fast protons. This kind of effect, which amounts to a rapid 
change of the heavy quark structure functions, can in principle be seen in 
high energy cosmic rays, as well as in LHC experiments. 

We would like to thank Carlos Pajares, Jaime Alvarez- Muniz, Nestor 
Armesto and J. Guilherme Milhano for discussions. This works has been 
done under contracts POCTI/36291/FIS/2000 and POCI/FIS/55759/2004 
(Portugal). 
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